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Simple  empirical  model  for  vibrational  spectra  of  single-wall  carbon  nanotubes 
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A  simple  empirical  model  and  approach  are  introduced  for  calculation  of  the  vibrational  spectra 
of  arbitrary  single  wall  carbon  nanotubes.  Differently  from  the  frequently  used  force  constants 
description,  the  model  employs  only  invariant  quantities  such  as  variations  of  lengths  and  angles. 
All  the  salient  qualitative  features  of  vibrational  spectra  of  nanotubes  naturally  follow  from  the 
vibrational  Hamiltonian  of  graphene  upon  its  isometric  mapping  onto  a  cylindrical  surface  and 
without  any  ad  hoc  corrections.  A  qualitative  difference  with  previous  results  is  found  in  a  parabolic, 
rather  than  a  linear,  long  wavelength  dispersion  of  the  transverse  acoustic  modes  of  the  nanotubes. 
The  parabolic  dispersion  is  confirmed  and  elucidated  in  the  provided  continuum  analysis  of  the 
vibrations.  We  also  discuss  and  use  an  alternative  definition  of  the  nanotube  unit  cell  with  only 
two  carbons  per  cell  that  illustrates  a  “true”  longitudinal  periodicity  of  the  nanotubes,  and  of  the 
corresponding  Brillouin  zone. 

PACS  numbers:  61.46.+W,  62.25.+g,  62.30.+d,  46.40.-f 


I.  INTRODUCTION 

Vibrational  spectra  of  individual  single  wall  carbon 
nanotubes  are  of  considerable  interest  and  have  been  cal¬ 
culated  previously  within  different  frameworks  such  as 
an  empirical  force  constant  model, 1-3  ab  initio  studies4 
and  tight-binding  molecular  dynamics.5  It  is  known  that 
the  higher-frequency  part  of  the  nanotube  spectra  is  rel¬ 
atively  well  represented  already  by  the  zone  folding  of 
the  graphene  spectrum.  The  lower-frequency  part,  how¬ 
ever,  has  generic  features  owing  to  the  one-dimensional 
character  of  nanotubes.  Particularly,  the  spectra  exhibit 
four  types  of  acoustic  modes  with  vanishing  frequencies: 
one  longitudinal,  two  transverse  and  one  twisting.  The 
existence  of  these  modes  has  to  do  with  general  consider¬ 
ations  -  displacements  of  the  tube  as  a  whole  along  and 
perpendicular  to  its  axis,  and  the  rotation  of  the  tube 
about  the  axis  do  not  cost  energy  -  rather  than  with  spe¬ 
cific  nanotube  interactions.  Analogous  vibrations  were 
also  discussed  in  the  context  of  quantum  wires  as  dilata- 
tional,  flexural  and  torsional  modes.6  These  modes  are 
important  contributors  to  the  low-temperature  quantized 
thermal  conductance  of  such  phonon  waveguides.7,8 

The  frequently  used  force  constant  model  of  Ref.  1  was 
first  developed  for  planar  graphene  based  on  the  experi¬ 
mental  data  for  graphite  and  then  adapted  for  nanotube 
geometries.  Direct  application  of  the  graphene  force  con¬ 
stant  values  was  not  found  to  lead  to  zero  frequencies 
for  all  four  modes  mentioned  above.  To  overcome  this 
difficulty,  special  curvature  corrections  were  introduced 
to  the  force  constants.1  In  this  paper  we  develop  another 
empirical  model,  where  the  harmonic  vibrational  Hamil¬ 
tonian  of  graphene  is  built  using  only  “invariant”  quan¬ 
tities  such  as  variations  of  bond  lengths,  interbond  and 
dihedral  angles.  Such  a  description  is  similar  in  spirit  to 
used  in  conformational  analysis  and  stereochemistry  (see, 
e.g.,  Ref.  9)  and  in  bond  models  for  vibrations  in  covalent 
semiconductors  (Ref.  10  and  references  therein).  Further 
mapping  of  the  graphene  Hamiltonian  onto  a  cylindrical 


surface  of  the  nanotubes  then  allows  to  derive  all  fea¬ 
tures  of  arbitrary  nanotube  spectra  “naturally” ,  without 
any  curvature  corrections.  This  way  the  idea  of  Ref.  1  of 
using  the  same  type  of  vibrational  Hamiltonian  for  both 
graphene  and  nanotubes  turns  out  to  be  realized  with  no 
need  for  ad  hoc  modifications. 

We  find  a  qualitative  difference  with  the  previously 
published  results  in  a  parabolic,  rather  than  a  linear,  dis¬ 
persion  of  the  transverse  acoustic  modes  of  nanotubes. 
This  parabolic  dispersion  is  further  illustrated  in  the 
corresponding  continuum  model  of  vibrations  and  is  in 
agreement  with  the  analysis  of  vibrations  of  elastic  cylin¬ 
drical  shells.11  The  continuum  model  also  shows  the  ori¬ 
gin  of  another  salient  long  wavelength  feature  of  the  nan¬ 
otube  spectra:  a  coupling  of  the  longitudinal  acoustic 
with  the  breathing  mode.  We  believe  the  parabolic  char¬ 
acter  of  the  low-frequency  part  of  the  transverse  mode 
dispersion  is  quite  generic  similarly  to  the  well-known 
bending  waves  of  rods12  with  wavelengths  much  longer 
than  the  rod  size.  Such  a  parabolic  dispersion  was, 
e.g.,  calculated  for  the  lowest  flexural  modes  of  quan¬ 
tum  wires;6,7  recent  applications  of  elastic  cylinder  mod¬ 
els  also  include  vibrations  of  cytoskeletal  filaments  and 
microtubules.13-15  To  our  knowledge,  the  parabolic  dis¬ 
persion  was  not  calculated  previously  for  carbon  nan¬ 
otubes. 

With  our  approach  we  can  easily  calculate  vibrational 
spectra  of  arbitrary  (TV,  M)  nanotubes.  In  doing  this,  we 
also  employ  an  alternative  definition  of  the  nanotube  unit 
cell  with  only  two  carbons  per  cell  as  opposed  to  possibly 
many  carbon  atoms  of  the  conventional  definition.1  This 
way  a  “true”  longitudinal  periodicity  of  the  nanotubes  is 
elucidated.  The  period  is  a  projection  of  one  of  the  prim¬ 
itive  vectors  onto  the  nanotube  axis,  the  primitive  vector 
itself  being  in  general  not  parallel  to  the  axis.  The  Bril¬ 
louin  zones  can  correspondingly  be  wider  than  usually 
used  and  the  total  number  of  the  vibrational  branches 
turns  out  to  be  just  6 dc,  where  dc  is  the  greatest  com¬ 
mon  divisor  of  N  and  M. 
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For  numerical  computations,  we  will  be  using  some 
parameterization  of  the  model,  specifically  based  on  a 
set  of  data  from  Ref.  1  as  well  as  on  some  experimental 
data.  However,  the  numerical  computations  here  serve 
mostly  illustrative  and  qualitative  purposes. 

II.  VIBRATIONAL  HAMILTONIAN 

Despite  the  fact  that  atoms  of  carbon  nanotubes  are 
arranged  in  a  3 -d  fashion,  the  excitations  of  nanotubes 
can  be  described  in  the  way  very  similar  to  excitations 
of  the  planar  graphene.  To  clearly  see  this  connection 
for  vibrational  excitations,  one  can  use  local  (position- 
dependent)  coordinate  systems  for  atomic  displacements 
and  the  elastic  potential  energy  written  in  an  “invariant” 
form.  This  way  the  qualitative  transformation  from  vi¬ 
brational  spectra  of  an  infinite  plane  to  spectra  of  curved 
cylindrical  structures  appear  naturally  without  ad  hoc 
corrections. 

Following  Ref.  1,  we  also  consider  carbon-carbon  elas¬ 
tic  interactions  up  to  the  fourth  nearest  neighbor.  The 
harmonic  potential  energy  U  =  U\  +  U0  is  however  ex¬ 
pressed  as  a  function  of  only  invariant  quantities  such  as 
variations  of  bond  lengths  and  various  angles.  The  first 
term  Ui  would  correspond  here  to  in-plane  deformations 
of  graphene  and  in  general  requires  ten  elastic  parameters 

K™- 

U  =  [Ki(Slij  +  5l?k)  +  K25(p?jk  +  K36hjSljk 
{ijk) 

+  +  51  jk)]  (la) 

4  {il}=m 

+  J2  £  [K^SkjSlki  +  K^S(pijkSipjki 

m=3  ( ijkl ) 

+  K™(5(pijk6lki  +  SipjkiSlij)] ,  (lb) 

while  term  U0  would  describe  out-of-plane  graphene  dis¬ 
tortions  and  needs  three  elastic  parameters: 

4  {iiy=m 

u0=J2  E  K™sx2m-  (2) 

m=2  {ijkl) 

The  structure  of  Eqs.  (1,2)  can  be  conveniently  thought 
of  in  terms  of  triangular  plaquettes  (ijk)  formed  by  bonds 
(ij)  and  (jk)  connecting  nearest  carbons  i  and  j,  and 
j  and  k,  respectively.  Variation  of  bond  (ij)  length 
is  denoted  Slij ,  and  variation  of  the  inter-bond  angle 
at  the  common  carbon  j  denoted  Stfijk-  Correspond¬ 
ingly,  Eq.  (la)  completely  describes  the  deformation  en¬ 
ergy  of  individual  plaquettes.  Equations  (lb, 2),  on  the 
other  hand,  completely  describe  the  interactions  of  de¬ 
formations  on  neighboring  plaquettes;  specifically,  (ijkl) 
stands  for  plaquettes  (ijk)  and  (jkl)  that  are  adjacent 
along  bond  (jk).  There  are  three  different  ways  to  form 
adjacent  plaquette  pairs  and  notation  {il}  =  m  distin¬ 
guishes  them  by  indicating  that  carbons  i  and  l  are  the 


FIG.  1:  Nomenclature  of  possible  triangular  plaquette  pairs 
adjacent  along  a  nearest-neighbor  carbon-carbon  bond.  To 
guide  the  eye,  individual  adjacent  plaquettes  are  filled  with 
lines  at  different  angles. 


mth  nearest  neighbors  (m  =  2,3,4),  as  illustrated  in 
Fig.  1.  (In  case  m  =  2,  it  is  actually  plaquettes  (ijk) 
and  (Ijk)  that  are  adjacent  along  (jk).  This  case  is  not 
explicitly  included  in  Eq.  (lb)  because  of  the  constraint 
that  a  sum  of  inter-bond  angles  for  three  plaquettes  sur¬ 
rounding  a  carbon  atom  is  fixed.)  The  “out-of-plane” 
interaction,  Eq.  (2),  involves  the  variations  dxijki  of  the 
dihedral  angles  between  corresponding  plaquettes. 

It  is  worth  stressing  that  for  the  planar  graphene, 
Eqs.  (1,2)  give  the  most  general  description  of  harmonic 
interactions  involving  up  to  the  forth  nearest  neighbors. 
As  such,  it,  of  course,  can  reproduce  a  force  constants 
description.  The  latter  would  be  derived  by  expanding 
Eqs.  (1,2)  in  the  atomic  pair  differences.  Evidently,  then 
there  would  be  certain  relationships  between  a  larger 
number  of  force  constants  as  dictated  by  the  invariance  of 
the  potential  energy  with  respect  to  overall  rotations.16 
Specifically,  ten  elastic  parameters  of  Eq.  (1)  yield  twelve 
force  constants,  comprising  eight  in-plane  constants  of 
the  type  explicitly  considered  in  Ref.  1  and  four  con¬ 
stants  mixing  radial  and  tangential  displacements  that 
were  implicitly  set  to  zero  in  that  reference.  Three  elas¬ 
tic  parameters  of  Eq.  (2)  yield  four  out-of-plane  force 
constants.1 

The  advantage  of  using  invariant  quantities  in 
Eqs.  (1,2)  is  that  the  same  functional  form  of  the  po¬ 
tential  energy  can  be  directly  used  when  carbon  atom 
positions  are  (isometrically)  mapped  from  the  graphene 
plane  onto  the  cylindrical  surface  of  a  nanotube.  When 
on  the  nanotube  surface,  variations  of  the  bond  lengths 
and  various  angles  just  need  to  be  calculated  from  the 
carbon  atom  displacements  using  the  actual  curved  ge¬ 
ometry.  Of  course,  a  new  set  of  force  constants  appro¬ 
priate  for  the  now  curved  geometry  can  be  again  derived 
through  the  pair  expansions.  Force  constants  so  obtained 
would  automatically  obey  the  correct  relationships  to  sat¬ 
isfy  the  invariance  with  respect  to  rotations. 

The  standard  translational  invariance  of  the  vibra¬ 
tional  Hamiltonian  is  preserved  if  displacements  of  car¬ 
bons  are  expressed  not  in  terms  of  common  (xyz)  co¬ 
ordinates  but  in  terms  of  local  orthogonal  coordinates 
( uvw ):  u  -  along  the  nanotube  axis,  v  -  perpendicular 
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to  the  tube  axis  and  parallel  to  the  tube  surface,  and  w 
-  perpendicular  to  the  tube  surface.  These  simple  ideas 
are  illustrated  in  more  detail  in  Appendix  A  for  an  easier 
to  follow  example  of  the  relationship  between  vibrations 
of  a  linear  chain  of  atoms  and  vibrations  of  a  ring  of 
atoms.  For  the  problem  at  hand,  we  use  the  local  dis¬ 
placements  (unai  vna,  wna)  depending  on  carbon  na  ac¬ 
tual  geometric  position  on  the  cylindrical  surface  of  the 
tube.  (Here  carbon  index  n a  consists  of  a  2 -d  vector  n 
specifying  the  unit  cell  of  the  parent  graphene  plane  and 
a  =  1,2  specifying  one  of  the  2  carbons  in  the  graphene 
unit  cell.)  With  local  displacement  bases  in  place,  the 
invariance  of  the  vibrational  Hamiltonian  with  respect 
to  translations  by  graphene  primitive  vectors  is  held  on 
equal  footing  in  graphene  and  nanotubes:  wave  vectors  k 
will  “know”  only  differences  between  neighboring  carbon 
indices  no.  The  problem  is  thereby  reduced  to  calcula¬ 
tion  of  the  usual,  “graphene-like”  ,6x6  dynamical  ma¬ 
trix,  but  which  naturally  contains  the  correct  mixing  of 
the  “in-plane”  and  “out-of-plane”  displacements  in  nan¬ 
otubes.  Using  proper  quantization  rules  for  the  allowed 
phonon  wavevectors,  one  can  then  readily  derive  the  vi¬ 
brational  spectra  of  nanotubes  of  arbitrary  chirality. 

III.  PRIMITIVE  CELLS  AND  BRILLOUIN 
ZONES 

The  geometry  of  a  single  wall  carbon  nanotube  is  de¬ 
termined  by  the  chiral  vector1 

Ch=Nai+Ma2,  (3) 

where  ai  and  a2  (|ai|  =  |  a2 1  =  a)  are  two  prim¬ 
itive  vectors  of  the  2-d  graphene  (hexagonal)  crystal 
structure.  Vector  Ch  is  perpendicular  to  the  nanotube 
axis.  Conventionally,1  the  translational  vector  parallel 
to  the  tube  axis  is  defined,  which  we  denote  here  as 
Td:  Td  =  t\ ai  +  <2&2,  where  t\  =  (2 M  +  7V)/cTr, 
f 2  =  —(27V  +  M)/du  and  du  is  the  greatest  common  di¬ 
visor  of  (27V  +  M)  and  ( 2M  +  7V).  The  resulting  unit  cell 
of  the  nanotube  built  of  Ch  and  To  can  contain  many 
carbons  7V0  =  4(7 V2  +  TV/2  +  NM)/dn,  and  the  longitu¬ 
dinal  period  T o  of  chiral  tubes  be  much  larger  than  a. 
The  corresponding  Brillouin  zones  (BZs)  would  then  be 
narrow  and  contain  many  excitation  spectrum  (whether 
vibrational  or  electronic)  branches. 

In  this  paper  we  use  an  alternative  picture  of  the  unit 
cell  and  BZ  construction  that  is  aimed  at  having  as  small 
number  of  branches  in  the  zone  as  possible.  As  described 
in  more  detail  in  Appendix  B,  this  number  of  branches  is 
determined  by  the  greatest  common  divisor  of  TV  and  TV/, 
denoted  by  dc.  The  total  number  of  continuous  branches 
in  BZ  is  equal  to  2 dc  per  each  degree  of  freedom  of  a  car¬ 
bon  atom,  that  is,  6 dc  for  vibrational  excitations.  This 
corresponds  to  the  nanotube  unit  cell  containing  only  two 
carbons  and  which  can,  e.g.,  be  built  with  primitive  vec¬ 
tors  Ch/dc  and  T.  Different  from  vector  T^,  the  trans¬ 
lational  vector  T  =  Pai  +  Qa-2  is  in  general  not  paral¬ 


lel  to  the  nanotube  axis.  As  discussed  in  Appendix  B 
(see  Eq.  (B7)),  integers  P  and  Q  here  satisfy  condition 
TV/P  —  NQ  =  dc.  It  is  the  projection  of  vector  T  onto 
the  axis  that  determines  the  longitudinal  period 

0|  |  =  V3dca2/2Ch  (4) 

and  the  width  2-k /a\\  of  the  BZ. 

The  number  of  branches  in  the  BZ  is  related  to  the 
transverse  quantization  of  the  2-d  wave  vector  k  of  the 
parent  graphene  band  excitations: 

kC  h  =  2t rZ,  (5) 

resulting  in  the  appearance  of  one-dinrensional  sub-bands 
characterized  by  the  integer  quantum  number  l.  The  con¬ 
struction  of  the  unit  cell  and  BZ  employed  in  this  paper 
recognizes  that  there  would  be  only  dc  “unique”  quan¬ 
tization  levels  (that  is,  only  dc  independent  integers  l ): 
all  other  allowed  k-vectors  can  be  obtained  with  trans¬ 
lations  by  graphene  reciprocal  vectors.  Many  branches 
of  the  conventional1  BZ  would  not  exhibit  gaps  at  that 
BZ  boundary;  they  would  correspondingly  become  sin¬ 
gle  continuous  bands  when  properly  “unfolded”  in  our 
construction.  The  quantization  lines  within  our  BZs  can 
span  several  hexagons  of  the  graphene  reciprocal  lattice 
(see  example  of  Fig.  9). 

IV.  VIBRATIONAL  SPECTRA 

As  was  mentioned  above,  potential  energy  in  Eqs.  (1,2) 
is  capable  of  reproducing  results  of  the  force  constant 
model1  for  the  graphene  spectrum.  Figure  2  (a),  how¬ 
ever,  has  been  calculated  with  a  parameterization  of  elas¬ 
tic  constants  A',  in  Eqs.  (1,2)  such  as  to  achieve  only  a 
close  similarity  to  the  published  spectrum.1  In  our  cal¬ 
culations,  we  chose  to  slightly  and  somewhat  arbitral- 
ily  modify  the  tangential  force  constants  (j>\  1  from  the 
published1  values  so  as  to  satisfy  +  6 +  Aej)'3'1  + 
14 =  0.  The  latter  equality  is  required  by  the  ro¬ 
tational  invariance  -  elastic  energy  should  be  zero  for 
the  overall  rotation  of  the  graphene  plane.  Original 
constants1  do  not  obey  it.  A  recently  published  new  set 
of  force  constants17  also  does  not  satisfy  this  require¬ 
ment.  The  overall  scaling  of  elastic  constants  was  chosen 
here  so  as  to  reproduce  graphene  experimental  optical 
frequencies  of  1580  and  868  cm-1.  Once  the  values  of 
elastic  constants  A',;  have  been  defined  for  graphene,  the 
same  values  are  used  to  calculate  spectra  of  nanotubes, 
examples  of  which  are  shown  in  Figures  2  (b)-(d).  As 
discussed  in  Sec.  II,  we  do  not  need  the  knowledge  of 
force  constants  of  Ref.  1  for  curved  geometries  because 
the  Hamiltonian  used  automatically  preserves  all  invari¬ 
ance  requirements. 

In  displaying  the  nanotube  spectra,  we  use  the  defi¬ 
nition  of  BZ  as  discussed  in  Sec.  Ill  and  Appendix  B. 
With  that  definition,  the  BZ  contains  only  6 dc  vibra¬ 
tional  branches.  For  armchair  tubes,  a||  =  a/2  and  the 
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FIG.  2:  Model  vibrational  spectra:  (a)  Graphene;  (b)  (10,10) 
tube;  (c)  (10,9)  tube;  (d)  (16,1)  tube.  See  text  for  definition 
of  aii . 


BZ  in  Figure  2  (b)  is  twice  as  wide  as  the  conventional 
BZ,  that  picture  would  be  restored  by  simple  folding. 
Figure  2  (b)  exhibits  dc  =  10  “independent”  transverse 
quantization  levels.  On  the  other  hand,  Figures  2  (c)  and 
(d)  with  dc  =  1  correspond  to  only  one  transverse  quan¬ 
tization  level,  that  is,  to  six  branches.  Each  polarization 
band/branch  exhibits  a  continuous  evolution  of  the  po¬ 
larization  vectors  with  in  general  a  strong  fc|  dependence. 
It  is  worth  noting  that  transverse  acoustic  modes  in  this 
picture  have  their  frequency  vanishing  at  finite  k\\  (which, 
of  course,  lie  in  a  hexagon  of  the  graphene  reciprocal  lat¬ 
tice  other  than  the  one  where  fcy  =  0  is). 

Figure  3  shows  low-frequency  parts  of  the  vibrational 
spectra  of  the  (10,10)  and  (10,0)  tubes  reduced  to  a  sin¬ 
gle  hexagon  of  the  reciprocal  lattice.  Panel  (a)  can  be 
directly  compared  to  the  published  results2  derived  from 
the  model  of  Ref.  1.  Apart  from  the  small  differences, 
such  as  values  of  acoustic  velocities,  likely  caused  by  our 
modification  of  the  force  constants,  there  is  one  qualita¬ 
tive  disparity.  The  dispersion  of  the  transverse  acoustic 
modes  in  Figure  3  is  clearly  seen  to  be  parabolic  in  con¬ 


FIG.  3:  The  low-frequency  part  of  calculated  nanotube  vibra¬ 
tional  spectra:  (a)  (10,10)  tube,  (b)  (10,0)  tube.  Note  that 
here  wave  vectors  fcy  are  measured  with  respect  to  closest  Y 
points  of  hexagons  of  the  reciprocal  graphene  lattice. 


trast  to  a  linear  dispersion  discussed  in  Refs.  1,2,5.  The 
parabolic  dispersion  is  seen  over  a  wider  range  of  fc||  on 
going  to  tubes  of  smaller  radii,  compare  Figures  3  (a) 
and  (b).  In  Sec.  V  we  give  an  analytic  confirmation  of 
this  observation  within  a  framework  of  a  continuum  me¬ 
chanics. 


V.  CONTINUUM  ANALYSIS 

Analysis  of  deformations  and  vibrations  of  thin-walled 
elastic  cylinders  goes  back  as  far  as  to  Rayleigh  and  Love; 
see,  e.g.,  Refs.  18-21,  references  therein,  and  Refs.  11,22 
for  a  dedicated  analysis  of  vibrations. 

In  the  case  of  a  planar,  graphene,  structure,  the  con¬ 
tinuum  potential  elastic  energy  can  be  written  as 

Ui  =  ^  j  dxdy[Ci{ux+vv)2 

+  C2  (( Uy  +  vx)2  -  4 uxVy))  ,  (6) 

UQ  —  —  J  dx  dy  [Di(wxx  -\-Wyy) 

+  D2{wly  -  WXXWyy)]  .  (7) 

Here  x  is  a  coordinate  that  would  later  become  along  the 
cylinder  axis  and  y  coordinate  along  the  cylinder  circum¬ 
ference,  p  is  the  mass  density.  Displacements  fields  u,  v 
and  w  would  become,  respectively,  parallel  to  the  cylin¬ 
der  axis,  parallel  to  its  circumference,  and  perpendicular 
to  the  cylindrical  surface;  x  and  y  subindices  denote  the 
differentiation  over  corresponding  coordinates.  This  type 
of  deformation  energy  is  well  known  in  the  continuum  me¬ 
chanics  of  plates12,20  and  can  be  readily  derived  from  the 
discrete  form  (1,2).  In  the  former  picture,  elastic  con¬ 
stants  Ci,  C2,  D\  and  D2  in  Eqs.  (6,7)  are  expressed  in 
terms  of  stretching  and  bending  rigidities  and  Poisson’s 
ratio.  In  the  latter  derivation,  they  would  be  expressed 
through  constants  IC s  in  Eqs.  (1,2).  Note,  however,  that 
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FIG.  4:  Branches  of  the  vibrational  spectra  obtained  in  the 
continuum  model  described  in  the  text  at  (a)  ky  =  0,  (b) 
ky  =  1  / R,  (c)  ky  =  n/R  with  n  >  2.  Calculations  were 
performed  for  R  =  6.78  A  that  one  would  have  for  the  (10,10) 
tube.  Elastic  parameters  used  in  calculations  for  this  Figure 
would  correspond  to  graphene  velocities  cia  =  24  km/s,  Cta  = 
14  km/s.  Dispersion  of  the  out-of-plane  graphene  vibrations 
u>  =  5k 2  was  taken  with  5  —  6  x  10-'  m2/s  (see  Ref.  4),  and 
D2 /D\  =  4C2  / Ci .  Graphene  results  are  shown  by  the  dashed 
lines  for  comparison.  Dispersion  of  the  graphene  transverse 
acoustic  mode  practically  coincides  with  that  of  the  nanotube 
twisting  mode. 


Eqs.  (1,2)  should  then  be  transformed  only  in  the  context 
of  purely  acoustic  deformations,  in  which  two  carbons 
in  a  hexagon  move  “in-phase”.  That  is  why  we  will  be 
dealing  here  with  the  3x3  matrix  in  Eq.  (9)  rather  than 
with  a  6  x  6  dynamical  matrix.  One  can  develop  a  con¬ 
tinuum  model  that  would  include  optical  “out-of-phase” 
deformations  as  well.  The  deformation  energy  of  type  of 
Eqs.  (6,7)  was  already  proved  to  be  useful  in  studies  of 
large  deformations  of  carbon  nanotubes.23 

The  modification  of  Eqs.  (6,7)  upon  formation  of  a 
cylindrical  body  in  continuum  mechanics  is,  e.g.,  dis¬ 
cussed  in  Refs.  11,20,21.  As  is  also  shown  in  Appendix 
A  for  a  discrete  model,  this  corresponds  to  a  simple  sub¬ 
stitution  in  (6,7):  vy  — >  vy  +  w/R  and  wy  — >  wy  —  v/R, 
R  being  the  cylinder  radius.  Using  this  substitution,  one 
can  easily  study  the  problem  of  small  vibrations  of  a  con¬ 
tinuum  cylinder.  The  vibrational  frequencies  w  for  the 
plane  waves  with  a  two-dimensional  wave  vector  (kx,ky) 
are  determined  from  the  eigenvalue  equation 

w2d  =  Md,  (8) 


where  displacement  vector  d  =  ( u,v,w )  and  matrix 
M(kx,  ky)  is 


M= 


'Cikl  +  C2k2y 
(Ci-C2)kxky 
— iAkx/R 


(Ci-C2)kxky  iAkx/R 

Ciky+C2k 2  —iBky/R 

iBky/R  uj2  +  Di(kl  +  ky)2_ 


(9) 

Here  C\  =  C\-\-D\j R2 ,  C2  —  C2-\-D2/ R2 ,  A  =  2C2  —  Ci, 
B  =  c1  +  {D1+  D2/2)k2x  +  Dik2  and  u2  =  C1/R2. 

In  the  case  of  the  planar  structure,  R  — >  00,  Eqs.  (8,9) 
lead  to  two  acoustic  waves  of  in-plane  vibrations  (with 


longitudinal  c;a  =  C\^2  and  transverse  c*a  =  C^2  veloc¬ 
ities)  and  an  acoustic  wave  of  the  out-of-plane  vibrations 
with  a  parabolic  spectrum  u>  =  D\^2(k2  +k2),  see  Figure 
4  (a). 

For  a  cylinder/tube  of  a  finite  radius  R ,  the  trans¬ 
verse  quantization  imposes  a  restriction  on  values  of 
ky  =  n/R ,  where  n  is  an  integer.  Of  special  interest 
to  us  here  are  values  of  ky  equal  to  0  and  to  ±1  / R. 
The  case  of  ky  =  0  (Fig.  4  (a))  yields  two  of  the  four 
acoustic  modes  of  the  tube  with  vanishing  frequencies: 
(i)  the  longitudinal  mode  with  a  low-frequency  disper¬ 
sion  coinciding  with  the  longitudinal  acoustic  wave  of 
the  graphene  and  (ii)  the  twisting  mode  whose  disper¬ 
sion  is  somewhat  modified  from  the  transverse  acous¬ 
tic  mode  of  graphene  by  virtue  of  the  elastic  constant 

~  i /2 

D2:  Cta  =  C2  .  The  longitudinal  mode  in  the  case 
of  the  cylinder  couples  with  the  breathing  mode,  whose 
frequency  at  kx  =  0  is  uib  and  whose  dispersion  is  deter¬ 
mined  by  constant  D\.  As  a  result  of  this  coupling,  an 
anti-crossing  behavior  of  the  branches  arises  as  is  clearly 
seen  in  the  Figure,  (iii)  The  two  other,  degenerate, 
acoustic  modes  of  the  tube  spectrum,  usually  referred 
to  as  transverse  acoustic  modes  for  carbon  nanotubes1 
or  as  flexural  modes  in  other  applications, 6,13,14,22  cor¬ 
respond  to  ky  =  ±1/R  (Fig.  4  (b)).  It  is  apparent 
from  the  Figure  that  these  modes  have  a  parabolic  spec¬ 
trum.  In  fact,  one  can  easily  show  this  analytically  by 
the  perturbation  analysis  of  (8,9)  in  kx:  contributions  to 
the  linear  coefficient  in  dispersion  oj(kx)  exactly  cancel. 
The  low-frequency  parabolic  dispersion  of  these  modes, 
although  modified  by  the  curvature  rigidity,  Eq.  (7), 
is  mainly  determined  by  the  in-plane  stretching  rigid¬ 
ity,  Eq.  (6).  Neglecting  Eq.  (7),  one  would  obtain  the 
long  wavelength,  for  kxR  <C  1,  dispersion  of  this  mode 
as  w(kx)  =  (2(Ci  —  C2)C2/Ci)112  k2R.  Evidently,  the 
parabolic  character  of  the  dispersion  becomes  even  more 
apparent  for  smaller-radius  tubes.  These  modes  corre¬ 
spond  to  bending  vibrations  of  the  tube  as  a  whole  and, 
in  this  sense,  are  similar  to  the  bending  vibrations  of 
rods,  whose  generic  long- wavelength  parabolic  dispersion 
is  well  known.12  Similarities  between  Figures  4  and  3  (a) 
are  obvious.  Note,  however,  some  quantitative  differences 
caused  by  different  values  of  effective  parameters. 


VI.  SUMMARY  AND  DISCUSSION 


We  have  described  a  simple  empirical  model,  in  which 
vibrations  of  the  graphene  and  individual  single  wall  car¬ 
bon  nanotubes  are  treated  and  calculated  on  the  same 
footing.  Differently  from  the  force  constant  model,  our 
model  uses  only  “invariant”  quantities:  variations  of 
bond  lengths,  interbond  and  dihedral  angles.  As  a  result, 
the  isometric  mapping  from  the  planar  graphene  onto  a 
cylindrical  surface  of  nanotubes  automatically  preserves 
all  the  right  relationships  between  equivalent  force  con¬ 
stants.  Importantly,  all  calculated  vibrational  spectra  of 
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nanotubes  correctly  exhibit  four  types  of  acoustic  exci¬ 
tations  with  vanishing  frequencies  (one  longitudinal,  two 
transverse  and  one  twisting) ,  which  are  derived  naturally 
and  follow  from  the  symmetries  of  the  underlying  system. 
Although  the  results  obtained  are  largely  similar  to  the 
earlier  published,  we  have  also  found  an  important  qual¬ 
itative  difference.  The  long-wavelength  dispersion  of  the 
transverse  acoustic  modes  is  shown  to  be  parabolic  rather 
than  linear.  One  consequence  of  this  is  that  the  vibra¬ 
tional  density  of  states  should  exhibit  a  one-dimensional 
singularity  near  zero  frequency.  We  cannot  exclude  that 
there  can  be  other  physical  implications  of  our  finding, 
although  this  apparently  is  not  the  case  for  the  quan¬ 
tized  ballistic  thermal  conductance,7,8  for  which  the  ex¬ 
act  dispersion  law  is  irrelevant.  These  long- wavelength 
(wavelength  much  larger  than  the  tube  radius)  transverse 
modes  correspond  to  bending  oscillations  of  a  nanotube 
as  a  whole  and,  therefore,  are  similar  to  bending  vibra¬ 
tions  of  rods,  whose  low-frequency  spectrum  is  known 
to  be  parabolic.12  The  parabolic  dispersion  is,  e.g.,  evi¬ 
dent  in  calculations  of  the  lowest  flexural  modes  of  quan¬ 
tum  wires.6,7  The  origin  of  the  parabolic  dependence  for 
nanotubes  has  been  analytically  illustrated  using  a  con¬ 
tinuum  elastic  model  similar  to  used  in  the  analysis  of 
vibrations  of  cylindrical  shells.11  The  continuum  model 
also  clarified  the  coupling  between  longitudinal  acoustic 
and  breathing  modes  resulting  in  the  anti-crossing  behav¬ 
ior  clearly  seen  in  the  calculated  spectra  of  nanotubes. 

The  simplicity  of  our  model  allows  us  to  easily  calcu¬ 
late  vibrational  spectra  of  nanotubes  of  arbitrary  chiral¬ 
ities.  To  better  handle  such  spectra,  we  employed  an 
alternative  definition  of  the  nanotube  unit  cell  with  only 
two  carbons  per  cell.  This  definition  reveals  a  “true” 
longitudinal  periodicity  of  carbon  nanotubes  that  can 
be  substantially  shorter  than  used  in  the  conventional 
definition.1 

Being  a  result  of  the  straightforward  isometric  map¬ 
ping,  the  model  presented  in  this  paper,  has  its  quali¬ 
tative  limitations.  For  instance,  it  does  not  automati¬ 
cally  yield  the  curvature- induced  softening  of  the  phonon 
modes  observed  in  ab  initio  studies  such  as  that  of 
the  twisting  mode,4  or  more  complex  dependencies  for 
the  breathing  mode.4,24  Of  course,  these  effects  are  not 
deducible  from  the  graphene  properties.  To  describe 
them,  the  model  would  need  to  include  both  the  relax¬ 
ation  of  nanotube  geometries  (mapping  would  not  be 
exactly  isometric)  as  well  as  an  explicit  dependence  of 
the  elastic  energy  on  nanotube  radius  and  with  the  axial 
anisotropy.27  We  discussed  such  a  generalization  in  the 
context  of  uniform  deformations,25  which,  e.g.,  leads  to 
the  chirality-dependent  stiffness  of  nanotubes.26  We  be¬ 
lieve  the  present  model  can  also  be  further  developed  in 
this  regard  and  correspondingly  parameterized  by  com¬ 
parison  with  ab  initio  calculations.  A  more  accurate  pa¬ 
rameterization  is,  needless  to  say,  required  even  for  the 
present  model.  Such  quantitative  aspects  have  not  been 
pursued  in  this  paper. 
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APPENDIX  A:  A  RING  EXAMPLE 

Here  we  study  in  more  detail  a  simpler  illustrative  ex¬ 
ample  of  the  relationship  of  vibrations  of  an  infinite  linear 
atomic  chain  and  a  ring  of  atoms.  For  clarity,  displace¬ 
ments  of  atoms  are  restricted  to  the  plane  in  which  a  ring 
and  chain  belong.  Atom-atom  interactions  correspond 
to  formation  of  bonds  and  result  both  in  stretching  and 
bending  rigidity  of  a  linear  chain.  Figure  5  is  a  picture 
of  the  ring  consisting  of  N  =  10  atoms,  which  is  assumed 
to  preserve  the  nature  of  interactions  in  the  chain.  The 
figure  shows  a  common  system  of  coordinates  (y,  z )  as 
well  as  local  systems  (v,  w )  related  to  the  tangential  and 
normal  displacements  of  the  corresponding  atom.  Obvi¬ 
ously,  the  kinetic  energy 

T  =  Y  +  *n)  =  y  (A1) 

n  n 

where  n  is  the  1  -d  positional  index  of  the  atom  along  the 
ring  circumference.  We  first  write  the  potential  energy 
of  a  linear  chain  in  an  invariant  form  as 

=  f  ZX  +  - (A2) 

n  n 

where  the  first  term  describes  the  stretching  and  second 
the  bending  rigidity.  Then  we  assume  that  the  same 
functional  form  holds  for  the  ring.  From  the  geometry  of 


FIG.  5:  A  ring  of  N  =  10  atoms  connected  by  interatomic 
“bonds”  (thicker  lines).  The  stretching  rigidity  corresponds 
to  variations  of  the  bond  lengths  Z,  the  bending  rigidity  to 
variations  of  the  inter-bond  angles  X-  Number  of  atoms  N 
(or  the  radius  of  the  ring)  determines  the  angle  ip  =  2n/N. 
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FIG.  6:  The  vibrational  spectra  of  the  linear  chain  (solid  lines) 
and  of  the  ring  (filled  circles)  of  atoms.  The  parameters  used 
for  the  plot  are  K/M  =  1  and  AT/M  =  0.1. 

Figure  5,  the  variation  of  the  bond  lengths,  as  a  function 
of  the  local  displacements,  is 

Sin  =  (vn+i-vn)  cos(ip/2)  +  (wn+1  +  wn)  sin(y>/2),  (A3) 
while  the  variation  of  the  inter-bond  angles 

l-Sxn  =  (2wn-wn+i-wn-i)cos(<p/2) 

+  (vn+1  -  vn-i)sm(ip/2).  (A4) 

For  the  linear  chain  (ip  =  0  and  vn  =  yn,  wn  =  zn), 
Eq.  (A3)  describes  a  conventional  longitudinal  stretching 
while  Eq.  (A4)  would  yield  a  conventional  local  curva¬ 
ture.  Tangential  and  normal  vibrations  in  a  linear  chain 
are  fully  decoupled.  The  corresponding  two  branches  of 
the  vibrational  spectrum  are  simply 

u,(k)  =  (4K/M)  sin2  (ka/2),  .  . 

w2(fc)  =  (16AVM)sin4(fca/2),  1  ’ 

where  a  =  l  is  the  distance  between  the  atoms  along  the 
chain.  This  spectrum  is  shown  in  Figure  6  with  solid 
lines. 

In  the  curved  system  (our  ring  with  p  ^  0),  on  the 
other  hand,  tangential  and  normal  displacements  are 
coupled,  as  is  well  known  in  the  elasticity  theory  for 
curved  surfaces12,20  and  clearly  seen  in  Eqs.  (A3,A4).  If 
we  were  to  use  the  common  system  of  coordinates  for  dis¬ 
placements  (yn,  zn),  the  atom  contributions  to  the  poten¬ 
tial  energy,  Eq.  (A2),  would  be  position  dependent  and 
the  translational  invariance  with  respect  to  n  — >  n  +  1 
would  be  lost.  With  the  local  coordinates  (vn,wn),  this 
invariance  is  preserved  and  one  can  directly  use  the  con¬ 
ventional  transition  to  k— states,  which  would  now  be 
quantized  as 

k  =  (2Tr/Na)i,  i  =  0,...N—l.  (A6) 

The  unit  cell  length  a  >  l  here  is  now  the  distance 
between  the  atoms  of  the  ring  along  its  circumference. 
The  derivation  of  the  spectrum  from  the  equations  of 
motion  follows  straightforwardly  from  Eqs.  (A1  A4), 


and  the  results  are  shown  in  Figure  6  with  filled 
circles.  In  accordance  with  Eq.  (A6),  it  is  now  a  set 
of  discrete  frequencies.  The  coupling  between  normal 
and  tangential  displacements  resulted  in  important 
qualitative  modifications  of  the  spectrum,  which  have 
been  obtained  exactly  and  naturally.  Particularly,  one 
notices  3  zero-frequency  modes.  Two  of  them  (with  i  =  1 
and  i  =  N  —  1,  or  i  =  —  1,  in  Eq.  (A6)  precisely  one 
wavelength  on  the  ring  circumference)  correspond  to  the 
displacements  of  the  ring  as  a  whole  in  two  orthogonal 
directions.  The  other  (with  i  =  0)  is  a  pure  tangential 
mode  describing  the  rotation  of  the  ring  as  a  whole.  The 
k  =  0  finite-frequency  mode,  on  the  other  hand,  is  a  pure 
normal,  breathing,  mode  that  “borrowed”  its  strength 
from  the  parent  stretching  oscillations  of  the  linear  chain. 

Equations  (A3,A4)  can  be  used  to  study  the  continuum 
limit,  when,  keeping  the  same  radius  R,  we  increase  the 
number  of  atoms  N  — >  oo  and  bond  length  l  — >  0.  Then, 
evidently, 

5l/l->[dv/dy  +  w/R],  , 

Sx/l  —  &  [dw/dy  —  v/R]  /dy.  '  ’ 

The  expressions  in  brackets  in  Eq.  (A7)  provide  a  recipe 
for  a  transition  from  the  continuum  model  of  a  chain  to 
the  continuum  model  of  a  ring.  The  subsequent  analysis 
of  vibrations  is  straightforward  with  the  wavevector  k 
quantized  as  kR  =  n  (n  being  an  integer)  yielding  three 
zero-frequency  modes  and  the  breathing  mode  as  in  the 
discrete  case  above. 

It  should  be  noted  here  that  in  reality  the  very  values 
of  the  unit  cell  length  a  and  elastic  constants  K ,  Kb  can 
in  fact  somewhat  differ  for  the  ring  and  the  linear  chain, 
and  the  difference  would  be  N— dependent  as  determined 
by  the  equilibrium  bond  length  in  the  ring.  In  this  sense, 
what  is  compared  in  Figure  6  is  the  structure  of  the  spec¬ 
tra  for  the  same  values  of  the  essential  parameters  upon 
isometric  mapping  of  the  linear  system  onto  a  circular 
ring.  In  addition,  the  elastic  energy  of  the  ring  could  in 
general  contain  terms  absent  in  the  energy  of  the  chain, 
such  as  SxnSln-  They  are  not  deducible  from  the  func¬ 
tional  form  of  the  linear  system  and  would  have  to  be 
explored  on  their  own.  We  also  note  that  Eq.  (A2)  can 
be  generalized  to  include  longer  range  interactions. 

APPENDIX  B:  TRANSVERSE  QUANTIZATION 
AND  UNIT  CELLS 

Two  carbons  connected  by  C h  (Eq.  (3))  on  the 
graphene  plane  correspond  to  the  same  carbon  on  the 
nanotube  after  wrapping.  The  effective  2 -d  cyclic  condi¬ 
tion  can  therefore  be  written  as 

P'.'; rn  T  C/j  —  Rnrn,  Rnm  —  TIS.I  T  ?77,<l2:  (Bl) 

Rnm  being  the  position  of  one  of  the  carbons  (there  are 
two  of  them)  of  an  arbitrary  unit  cell  of  graphene.  If  k 


FIG.  7:  Unwrapped  (3,2)  tube  with  carbons  rearranged  to 
form  a  horizontal  strip  here.  The  tube  axis  is  shown  by  the 
dash-dotted  line.  Here  dc  =  1  and  the  T  vector,  Eq.  (B2), 
is  defined  by  the  ( P,Q )  pair  (2,1).  Small  circles  connected 
by  thinner  lines  indicate  a  sequence  of  carbons  reached  by 
the  T— translation.  It  is  transparent  that  all  hexagons  of  the 
strip  would  be  visited  this  way,  forming  an  effectively  1-d 
enumeration  of  Eq.  (B5).  The  conventional  vector  Td  = 
7ai  —  8a2  and  such  a  unit  cell  would  have  38  hexagons. 

is  the  2 -d  wave  vector  of  the  band  excitations,  then  in 
the  infinite  graphene  plane  it  would  have  two  indepen¬ 
dent  continuous  components.  In  nanotubes,  the  cyclic 
condition  leads  to  the  transverse  quantization  of  Eq.  (5). 
In  other  words,  wave  vectors  k  =  +  k||  allowed  by 

Eq.  (5)  lie  only  on  certain  quantization  lines  in  the  recip¬ 
rocal  plane  of  graphene,  kx  being  quantized  according 
to  Eq.  (5)  and  k|j  being  actually  a  1-d  continuous  wave 
vector  parallel  to  the  tube  axis. 

Let  us  define  vector 

T  =  Pan  +  Q&2  (B2) 

with  integer  P  and  Q.  From  Eqs.  (3)  and  (B2),  one  finds 
ai  =  (~QCh  +  MT)/ A,  a2  =  (PCh  -  NT) /A,  (B3) 
where 

A  =  MP  -  NQ.  (B4) 

Evidently,  if  one  can  find  such  P  and  Q  (of  course,  we  are 
interested  in  the  “smallest”  P  and  Q)  for  a  given  tube 
( N,M )  that  A  =  ±1  in  Eq.  (B4),  then  original  primi¬ 
tive  vectors  ai  and  a2  in  Eq.  (B3)  will  be  represented 
through  integer  amounts  of  C h  and  T.  Correspondingly, 
an  arbitrary  vector  R„m,  Eq.  (Bl),  will  be  expressed 
through  integer  quantities  of  C h  and  T  as  well.  Since  the 
cyclic  condition  defines  Rnm  in  Eq.  (Bl)  with  accuracy 
to  C h,  this  would  actually  mean  that  the  position  of  any 
hexagon  in  the  unwrapped  tube  is  determined  through  a 
single  vector  T: 

R„m  — >  Rj  =  iT,  (B5) 

where,  e.g.,  i  =  nM  —  mN  for  A  =  1.  One  can  think  of 
the  corresponding  unit  cell  built  of  and  T  that  would 
contain  only  2  carbons  in  the  cell.  Of  course,  vector 
T  does  not  have  to  be  parallel  to  the  tube  axis.  The 


FIG.  8:  Unwrapped  (4,2)  tube  with  carbons  rearranged  to 
form  a  horizontal  strip  here.  The  tube  axis  is  shown  by  the 
dasli-dotted  line.  Here  dc  =  2  and  the  T  vector,  Eq.  (B2),  is 
defined  by  the  “reduced”  ( P,Q )  pair  (1,0).  Two  sequences  of 
carbons  resulting  from  T— translations  are  shown  by  thinner 
lines:  small  circles  of  one  sequence  are  connected  by  the  solid 
line,  and  circles  of  the  other  by  the  dashed  line.  Translation 
from  one  sequence  to  the  other  is  achieved  with  vector  Ch/2. 
All  hexagons  of  the  strip  would  be  visited  this  way,  forming 
effective  enumeration  of  Eq.  (B8).  The  conventional  vector 
Td  =  4ai  —  5a2  and  such  a  unit  cell  would  have  28  hexagons. 

projection  of  T  on  the  nanotube  axis  7j  =  \/3a2  A/2Ch 
-  it  would  be  directed  in  opposite  ways  for  A  =  1  vs 
A  =  —  1  and  its  modulus  determines  the  corresponding 
longitudinal  period 

a||  =  V3a2/2Ch ,  (B6) 

where  Ch  =  |C/j|.  Note  that  vector  T  is  different  from 
the  symmetry  vector  defined  in  Ref.  1. 

It  is  easy  to  see  that  the  picture  described  in  the  pre¬ 
vious  paragraph  is  indeed  realized  whenever  dc  =  1  -  we 
will  call  it  the  irreducible  case  -  where  dc  is  the  greatest 
common  divisor  of  N  and  M.  (Considering  only  posi¬ 
tive  N  and  M  does  not  restrict  the  generality.)  Then, 
in  fact,  any  integer  value  of  A  in  Eq.  (B4)  can  be  estab¬ 
lished  with  an  appropriate  choice  of  P  and  Q.  Examples 
of  positive  (P,  Q)  pairs  satisfying  A  =  1  are  listed  here 
as  (N,M)  -  (P,Q):  (3,1)  (1,0),  (3,2)  -  (2,1), 

(5,2)  — >  (3,1),  and  (10,9)  — >  (9,8).  An  illustration  for 
the  (3,2)  tube  is  shown  in  Figure  7.  A  visual  picture  of 
the  irreducible  case  is  that  of  a  1-d  chain  (of  period  |T|) 
that  is  wrapped  around  the  nanotube  cylinder  with  an 
appropriate  helix  angle,  as  is  clearly  seen  from  Figure  7. 

“Reducible”  cases  -  with  dc>  1  -  can  be  described  in  a 
similar  but  somewhat  more  involved  fashion.  One  would 
extract  the  irreducible  structure  factors  Ni  and  Mp 

N  =  dcN1,  M  =  dcMi, 

so  that  the  greatest  common  divisor  of  Ni  and  Mi  equals 
1.  Then  the  procedure  described  above  for  the  irreducible 
case  can  be  applied  for  the  geometry  (iVi,Mi)  resulting 
in  the  first  translational  vector  T.  In  other  words,  inte¬ 
gers  P  and  Q  of  Eq.  (B2)  should  satisfy  condition 

MP  -  NQ  =  dc.  (B7) 

One  however  could  not  visit  all  hexagons  by  using  only 
so  defined  T.  The  needed  second  translational  vector 


9 


FIG.  9:  An  illustration  of  the  range  of  variation  of  k||  for 
(3,2)  and  (4,1)  tubes.  In  the  former  case  K2  =  2bi  —  3b2, 
in  the  latter  case  K2  =  bi  —  4b2 .  The  corresponding  ranges 
are  shown  as  thick  lines  between  the  centers  of  hexagons  of 
the  reciprocal  lattice  of  graphene.  Each  line  crosses  several 
hexagons.  If  all  inside-a-hexagon  segments  are  displaced  to 
one  hexagon  with  a  graphene  reciprocal  vector,  they  would 
form  a  traditional  picture  of  quantization  lines  in  the  first  BZ 
of  graphene. 

can  be  found  as  C h/dc.  The  resulting  enumeration  of  all 
hexagons  will  then  read  as 

R„m  ->  Ry  =  iT  +  jCh/dc,  j  =  0, 1, . . .  dc  -  1.  (B8) 

Index  i  here  is,  as  in  Eq.  (B5),  an  arbitrary  integer 
that  would  define  unique  carbons  while  j  results  only 
in  dc  unique  translations  due  to  the  cyclic  condition  in 
Eq.  (Bl).  An  illustration  for  the  (4,2)  tube  is  shown  in 
Figure  8.  A  visual  picture  of  the  irreducible  case  is  then 
of  dc  1  -d  chains  (of  period  |T|)  that  are  wrapped  around 
the  nanotube  cylinder  with  an  appropriate  helix  angle. 
The  elementary  unit  cell  contains  two  carbons  and  built 
of  vectors  T  and  C h/dc.  Longitudinal  periodicity  will 
again  be  determined  by  the  projection  of  vector  T  on 
the  nanotube  axis  which  is  Eq.  (4)  becoming  Eq.  (B6)  at 
dc  =  1. 

For  the  reciprocal  lattice  vectors,  defined  through  Ki  • 
C h  =  27t,  K2  ■  T  =  27T,  Ki  •  T  =  0,  and  K2  •  C h  =  0,  one 
readily  obtains 

=  —Qhi  +  Pb2,  K2  =  M1b1-JV1b2) 

where  bi  and  b2  are  primitive  graphene  reciprocal  vec¬ 
tors  defined  with  respect  to  ai  and  a2.  One  easily  finds 
that  |  K2 1  indeed  coincides  with  27r/a||.  An  illustration 
for  cases  of  (3,2)  and  (4,1)  tubes  is  given  in  Figure  9. 

This  way  one  arrives  at  the  band  excitation  picture 
with  2 dc  (for  each  component  of  the  wave  function  or 
polarization)  continuous  bands  corresponding  to  2  atoms 
in  the  unit  cell  and  dc  transverse  quantization  levels,  e.g. 
with  quantum  numbers 

l  =  —dc/2  +  1, . . . ,  0, . . .  dc/2.  (B9) 


The  latter  define  the  quantized  component  of  the  parent 
2-d  quasi-momentum  k  perpendicular  to  the  tube  axis: 

k±  =  (2tt /Ch)l.  (BIO) 

The  continuous  component  k\\  parallel  to  the  tube  axis, 
on  the  other  hand,  would  be  defined  within  a  BZ  whose 
width  is 

27r/a||  =  4nCh/VSdca2 .  (Bll) 

One  can  easily  calculate  that  the  resulted  number  of 
continuous  bands  and  the  width  of  the  BZ  lead  to  the 
correct  total  number  of  states  in  the  system,  which 
corresponds  to  2Ch/V3a2  hexagons  per  unit  length  of 
the  nanotube. 

The  conclusion  that  there  are  only  dc  “truly  unique” 
quantization  levels  can  be  confirmed  another  way  as  well. 
If  G  =  nibi  +n2b2  is  the  vector  of  the  reciprocal  lattice 
of  graphene,  then  ki  =  k  +  G  is  physically  equivalent 
to  k.  In  general,  k  and  ki  may  correspond  to  different 
quantum  numbers  in  Eq.  (5),  say  l  and  l\.  One  derives 
the  integer  difference  of  these  quantum  numbers  as 

81  =  l  —  l±  =  n\N  +  n2M  =  dc(n\Ni  +  n2TWi).  (B12) 

Once  again,  the  last,  “irreducible”,  factor  in  Eq.  (B12) 
can  take  any  integer  values: 

nilVi  +  n2Mi  =  0,  ±1,  ±2, . . . , 

with  appropriate  choices  for  integer  ni  and  n2.  It  is  then 
clear  that  the  number  of  unique  sub-bands  in  the  ex¬ 
tended  scheme  equals  precisely  dc  -  one  can,  e.g.,  choose 
quantum  numbers  of  Eq.  (B9)  for  “independent”  sub¬ 
band  indexing,  all  other  values  of  l  would  be  reducible 
to  independent  values  with  the  appropriate  adjustment 
for  the  one-dimensional  continuous  quasi-momentum  fc|| . 

It  is  worth  mentioning  that  in  the  ( TV,  M )  family  of 
nanotubes  with  TV  >  M ,  it  is  the  armchair  and  zigzag 
nanotubes  that  have  maximal  number  of  transverse 
quantization  levels  dc  =  TV.  The  tubes  immediately 
next  to  them,  (TV,  TV  —  1)  and  (TV,  1),  on  the  other  hand, 
would  have  only  one  quantization  level.  The  armchair 
tubes  have  then  the  longitudinal  period  of  a/2,  twice  as 
small  as  is  with  the  conventional  definition.1  One  can 
easily  see  it  with  T  vector  defined  by  the  (P,  Q)  pair  (1,0). 
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